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Abstract 

We present a supersymmetric SO(10) grand unified theory (GUT) of flavor based on an S4 family 
symmetry. It makes use of our recent proposal to use SO (10) with type II seesaw mechanism for 
neutrino masses combined with a simple ansatz that the dominant Yukawa matrix (the 10-Higgs 
coupling to matter) has rank one. In this paper, we show how the rank one model can arise 
within some plausible assumptions as an effective field theory from vectorlike 16 dimensional 
matter fields with masses above the GUT scale. In order to obtain the desired fermion flavor 
texture we use S4 flavon multiplets which acquire vevs in the ground state of the theory. By 
supplementing the S4 theory with an additional discrete symmetry, we find that the flavon vacuum 
field alignments take a discrete set of values provided some of the higher dimensional couplings 
are small. Choosing a particular set of these vacuum alignments appears to lead to an unified 
understanding of observed quarkdepton flavor: (i) the lepton mixing matrix that is dominantly 
tri-bi-maximal with small corrections related to quark mixings; (ii) quark lepton mass relations 
at GUT scale: m& — m T and m^ ~ 3m s and (iii) the solar to atmospheric neutrino mass ratio 
"JeMatm — #Cabibbo in agreement with observations. The model predicts the neutrino mixing 
parameter, U e s ~ #Cabibbo/ (3v2) ~ 0.05, which should be observable in planned long baseline 
experiments. 
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I. INTRODUCTION 



A unified understanding of the diverse pattern of quark lepton masses and mixings is a 
fundamental challenge for physics beyond the standard model l|. The two major elements 
of this flavor puzzle that any theory must explain are: (i) strong mass hierarchy in the 
quark and charged lepton sector and weak hierarchy for neutrinos; (ii) large lepton mixings 
i.e. 6*23 ~ 45° and 9[ 2 — 35° as against small quark mixings 9l 3 ~ 2.5° and 6\ 2 ~ 13° and 
apparent relation between some of the mixing angles and the fermion masses. Since grand 
unified theories (GUT) not only unify different gauge couplings at a high scale but also 
unify quarks and leptons within a single framework, they have often been thought of as an 
attractive venue for unraveling this puzzle. Furthermore the fact that the seesaw mechanism 
for understanding small neutrino masses Q also seems to require a B — L breaking scale close 
to the scale of coupling unification, makes this suggestion quite promising. The constraints 
of higher symmetry however make it highly nontrivial to understand all the details of flavor 
puzzle although many attempts have been made [sj]. 

In a recent paper, we have suggested a possible way Q to address this problem in su- 
persymmetric SO (10) GUT models. The main assumptions of ref. ^| are: (a) all fermion 
masses arise from effective Yukawa couplings [ji] involving 10 and 126 Higgs multiplets; (b) 
neutrino masses arise (J from type II seesaw mechanism and (c) the 10-Higgs Yukawa 
dominates fermion masses and has rank one. We showed in ref. 4} how this program when 
implemented using the already mentioned Higgs content of a single 10, 126 plus possibly 
another 10 or 120 Higgs fields not only explains all the qualitative features of quark and 
lepton flavor noted above but also makes a prediction for the lepton mixing angle U e % or # 13 . 
In most models we discussed in the apparent tri-bi-maximal mixing pattern observed 
for neutrinos did not arise from any symmetry. In this note, we pursue program outlined 
in j^J further by using discrete family symmetries to make this ansatz more predictive. Our 
strategy is to use flavon fields whose vevs give the effective Yukawa couplings responsible for 
fermion masses at the GUT scale. We use additional discrete family symmetries whose role 
is to constrain the ground state of the flavon Hamiltonian such that they lead to particular 
textures for the fermion mass matrices within certain assumptions. We are able to isolate 
a set of allowed flavon vacuum states which are such that the dominant part of the lepton 
mixing matrix naturally has a tri-bi-maximal form, provided some of the higher dimensional 
terms in the flavon superpotential are small. The desired flavon vacuum alignment seems 
to arise naturally with an S4 symmetry jol which unifies all three families of fermions into a 
32 multiplet. 

The new results of this paper are : (i) we show how the rank one model can arise naturally 
as an effective field theory from vectorlike 16 dimensional matter fields with masses above 
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the GUT scale; and (ii) how the detailed fermion flavor textures arise from the vacuum field 
alignments of gauge singlet S4 flavon fields leading to the following results naturally without 
adjustment of parameters: (a) the lepton mixing matrix has dominantly tri-bi-maximal form 
with small corrections related to quark mixings; (b) quark lepton mass relations at GUT 
scale: mb ~ m T and m M ~ 3m s and (c) the solar to atmospheric mass ratio m Q /m atm ~ 
#Cabibbo m agreement with observations. 



II. OVERVIEW OF THE SUSY SO(10) RANK ONE STRATEGY 

We use the Higgs fields that give fermion masses to consist of two 10 dimensional mul- 



tiplets (denoted by H,H') and a single 126 + 126 (denoted by A and A). The Yukawa 
superpotential for this case in a generic SO(10) model can be written as: 

W Y = hipi/jH + f ipipA + ti ipi/j H' , (1) 

where the symbol ip stands for the 16 dimensional representation of SO(10) that represents 
the matter fields. The coupling matrices h, h! and / are symmetric. As we show later in this 
paper, their detailed texture will be determined by the S4 symmetry. The representations 
H, H' and A each have two standard model (SM) doublets in them. The general way to 
understand how the two MSSM doublets arise from them is as follows: at the GUT scale 
Mj/, after the GUT and the B — L symmetries are broken, one linear combination of the 
up-type doublets and one of down-type ones remain almost massless whereas the remaining 
ones acquire GUT scale masses just like the color triplet and other non-MSSM multiplets. 
The electroweak symmetry is broken after the light MSSM doublets (to be called H u ^) 
acquire vacuum expectation values (vevs) and they then generate the fermion masses. The 
resulting formulae for different fermion masses are given by: 

Y u = h + r 2 f + r 3 ti, (2) 

Y d = ri (h + f + h'), 

Y e = ri (h-3f + c e h'), 
Y u d = h 3r 2 / + c u ti, 

where Y a are mass matrices divided by the electro-weak vev v w k and ri and c e>l/ are the 
mixing parameters which relate the H u d to the doublets in the various GUT multiplets. 
More precisely, the matrices h, f and h! in Y a are multiplied by the Higgs mixing parameters 
when they appear in the fermion mass matrices. The definitions of the couplings and the 



Higgs mixing parameters are given in ref. 1QJ. In our particular case with a second 10 



Higgs (H'), c e = 1 and c„ = r 3 . Furthermore, we use the type II seesaw formula for getting 
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neutrino masses which is possible to obtain with symmetry breaking pattern in SO(10) as 
given in 11]. 



M v = fv L . (3) 

Note that / is the same coupling matrix that appears in the charged fermion masses in Eq. 
(j2J), up to factors from the Higgs mixings and the Clebsch-Gordan coefficients. This helps 
us to connect the neutrino parameters to the quark-sector parameters. The equations (J2]) 
and (|3J) are the key equations in our unified approach to addressing the flavor problem. 
The main hypothesis of our approach in ref. {^J is that 

• the fermion mass formula of Eq. (j2J) are dominated by the matrix h with the contri- 
butions of / and hi being small perturbations; 

• the matrix h has rank one. 

It follows from these assumptions that in the limit of /, hi — > 0, the quark and lepton 
mixings vanish as do the neutrino masses. Once /, hi are turned on, one can choose / to 
be diagonal by an appropriate choice of basis and without any loss of generality. Since the 
neutrino masses are diagonal in this basis, the entire Pontecorvo-Maki-Nakagawa-Sakata 
(PMNS) matrix comes from the matrix that diagonalizes the charged lepton mass matrix 
and for arbitrary form of the later, the PMNS matrix will in general have large mixing 
angles. On the other hand, the Cabibbo-Kobayashi-Maskawa (CKM) matrix Vckm — U^U^ 
which in the limit of /, hi — Y is equal to a unit matrix, owes the origin of quark mixings 
to f,h'. The quark mixings are proportional to and hence small as observed. It is 

also clear that the charged lepton and quark masses of second and first generation are also 
proportional to and thus hierarchical. 

Our procedure in this paper is as follows: we supplement the above rank one hypothesis 
by a discrete family symmetry S4 so that forms of the h, f, hi are consequences of the 
vacuum expectation values of gauge singlet but S4 non-singlet flavon fields, 0$, thereby 
making the model more predictive. To implement this procedure, we first derive the GUT 
scale effective Lagrangian from a pre-GUT scale theory that has vectorlike 16-dim. matter 
spinor with masses slightly above the GUT scale. The resulting effective theory involves 
non-renormalizable higher dimensional operators involving ip, Higgs fields and the flavon 
fields <pi whose vevs generate the flavor texture observed at GUT scale. These are then 
extrapolated to the weak scale to compare with observations. 
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TABLE I: The fields and representations to generate the desired Yukawa couplings, uj = e l ~ . 



III. S 4 FAMILY SYMMETRY AND MODEL OF FLAVOR 

The S4 group is a 24 element group describing permutations of four distinct objects and 
has five irreducible representations with dimensions 3i©3 2 ©2©l 2 ffili. The distinction 
between the representations with subscripts 1 and 2 is that the later change sign under the 
transformation of group elements involving the odd number of permutations of S4. For other 
details of S4 group, see js]. 

We assign the three families of 16-dim. matter fermions ip to 3 2 -dim. representation of 
54 and the Higgs field H, A and H' to li, 1 2 , and li reps, respectively. We then choose 
three SO(10) singlet flavons 0» transforming as 3 2 ,3i,3 2 reps of S4 and one gauge and ^4 
singlet fields Si,S2 transforming as 1 2 and li respectively. We further assume that at a 
scale slightly above the GUT scale, there are two S4 singlet vectorlike pairs of 16 © 16 fields 
denoted by ipy and ipy. In order to get the desired Yukawa couplings naturally from this 
high scale theory, we supplement the 64 group by an Z n group with all the above fields 
belonging to representations given in the Table I. 

The most general high scale Yukawa superpotential involving matter fields invariant under 
this symmetry is given by: 

W = (0iV>)t/Vi + 4>vii>viH + Mr^rVvi (4) 
1V1 p 

+ j^s 2 ((f> 3 ^)A + —(faW)H', 

where the brackets stand for the S4 singlet contraction of flavor index. The singlet field Sj 
can have large vev as follows: consider its Z n charge to be such that the only polynomial 
term involving the Sj in the superpotential has the form s^/Mp -3 (in order to describe the 
essential potential, we ignore a possible s^s^ 2 term). The dominant part of the potential in 
the presence of SUSY breaking has the form: 

V( Si ) = -ml\ Sl \ 2 + k^ ZE + ---. (5) 
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Minimizing this leads to (sj) ~ [m 




i -4 , which is above GUT scale for larger values 



of the integer ki (which in turn is determined by the Z n symmetry charge of Sj). One could 
also have large vevs for s±, S2 by using anomalous U(l) charges for them using D-terms to 
break the U(l) symmetry. 

The effective theory below the scales M± t 2 and (sj) of the vector-like pair masses and the 
Sj-vevs respectively is given by: 



where we have omitted the dimensional coupling constants to make it simple for the purpose 
of writing. The discrete symmetries prevent 2 /M 2 corrections to these terms. So our pre- 
dictions based on this effective superpotential do not receive large corrections. We note that 
the non-renormalizable terms in Eq.(j4]) can also be obtained from renormalizable couplings 
if we introduce further ^-triplet vectorlike fields. Here, however we use only S^-singlet 
vectorlike fields to get rank 1 contribution to h and / Yukawa couplings and that is why 
we need the non-renormalizable terms to be present in Eq.(jH A few comments are in order 
regarding the need for the extra Z n symmetry. 

• The Z n group provides a selection rule of the flavon couplings and the charges of various 
fields under this are chosen so as to forbid direct renormalizable Yukawa coupling, e.g., 
(ijjijj)H, which can lead to loss of rank one property and hence the hierarchy of fermion 
masses. 

• The barred flavon fields 0j are introduced to obtain the potential of the flavons neces- 
sary for our vacuum alignment. They do not couple to matter fields. 

• We note that the replacement of 0i with 03, 03 is forbidden if a — b ^ 0, —4, and 
similarly unwanted terms can be forbidden when n is a large number. 

• The term 1 -0'0AS'i is Z n invariant, but transforms as I2 under S4 because ipip is 
symmetric due to SO(10) algebra and thus it is not allowed either. 

• The 1S4 invariant singlet S2 is introduced to forbid 0203, 0203 terms, which are un- 
wanted in the flavon superpotential. 

• The Z n symmetry allows mixed higher dimensional terms of the form 0j0j0j0j terms 
with i 7^ j. We assume that the couplings of these terms are small compared to other 
terms so that the alignment shift caused by these terms compared to that given below 
is small and does not affect our result. 



W = {<j>^){<j>^)H + (0 2 V) (0 2 VO A + (03^) A + {<fcH)H', 



(6) 
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The details of the flavon superpotential will be discussed later. 

In order to get fermion masses, we have to find the alignment 12j of the vevs of the flavon 
fields 01,2,3- We show below that the following choice of vevs are among the minima of the 
flavon superpotential provided the couplings of mixed terms between different 0j's are small 
compared to other couplings: 



V 1 / 
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1 
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Clearly, there are other vacua for the flavon model that we do not choose. What is however 
nontrivial is that the alignments are along quantized directions. This is a consequence of 
supersymmetry combined with discrete symmetries in the theory. Given these vev, we find 
from Eq. ([6]) that the Yukawa coupling matrices h, f, h' have the form: 
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(10) 



and the charged fermion mass matrices can then be inferred. The neutrino mass matrix in 
this basis has the form: 
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(11) 



/ 



where c/a = A C 1. It is diagonalized by the tri-bi-maximal matrix 

( ^ fi o ^ 

Utb = 



\ 



(12) 



/ 



This is however not the full PMNS matrix which will receive small corrections from diag- 
onalization of the charged lepton matrix, which not only make small contributions to the 
8 atm and #0 but also generate a small #13. 



The neutrino masses are given by m 



v3 



2a — c ; m v2 = 2c and m v \ 



c. To fit 



,2 

"a tin 



0.2, which fixes the neutrino masses 



observations, we require A = c/a ~ \J Am / Ar 
m„3 ~ 0.05 eV, m v2 ~ 0.01 eV, and m v \ ~ 0.005 eV. We will see below that A is also the 
Cabibbo angle substantiating our claim that neutrino mass ratio and Cabibbo angle are 
related. 

For the charged lepton, up and down quark mass matrices, we have: 
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(13) 



M, 



tan (3 
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M„ 
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nil + ^ mo —m® + mi 
nil — 5 —Trio + m i m o + M 

r 2 nii + r 3 5 r 2 nii — r 3 5 
r 2 m + r 3 nii 
r 2 niQ + M 





r 2 nii + r 3 5 r 2 m 
r 2 nii + r 3 5 -r 2 m + r 3 nii 



where tan /3 is a ratio of H u d vevs. Note that nii/niQ = A ~ 0.2 and of course niQ <C Af. A 
quick examination of these mass matrices leads to several immediate conclusions: 



1. The model predicts that at GUT scale m;, 

2. Since (M d ) n -> 0, we get V us ~ ^/m d /i 



m. 



3. The empirically satisfied relation m^m e ~ m s m d can be obtained by the choice of 
parameters — 3mi + 5 = (nil + 5)e lcr , where a is a phase. Solving this equation, we 
find that 5 = nii(l + icotcr/2). We obtain V us ~ (1 — r 3 /r 2 )5/m , thereby relating 
Cabibbo angle to the neutrino mass ratio m /m atm ~ A. 



4. m M ~ 



-3m 9 . 



5. The leptonic mixing angle to diagonalize Mi is related to quark mixing 9[ 2 ~ \V US , 



which leads to a prediction for sin 13 



U e3 ~ ¥*. ~ 0.05 



3^ 



13|- 



6. Kb ~ — cot 9 



atm- 
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7. The masses of up and charm quarks are given by the parameters ^3 and are therefore 
not predictions of the model. 

8. CP violation in quark sector can put in by making the parameters h' complex. 

9. The model predicts a small amplitude for neutrino- less double beta decay from light 
neutrino mass: m Vee ~ csin 6[ 2 ~ 0.3 meV. 

The first four relations are fairly well satisfied by observations; the fifth prediction (i.e. 
that for Uez) can be tested in upcoming reactor and long baseline experiments. Note that 
the deviation from tri-bi-maximal mixing pattern coming from the charged lepton mass 
diagonalization could be thought of as a small perturbation of the neutrino mass matrix [3] 
except that we predict the form of the perturbation from symmetry considerations. The 
sixth prediction gives a smaller value for V c b (0.02 as against observed GUT scale value of 
0.03) if one uses GUT scale extrapolated value of the known b mass. However, in the MSSM 
there are threshold corrections to the b— s quark mass mixing from gluino and wino exchange 
one-loop diagrams; by choosing this contribution, one could obtain the desired 

Note that in this model, the top quark Yukawa coupling at GUT scale arises from an 
effective higher dimensional operator. We have showed the effective operator in Eq.(j6]) by 
expanding <j)/M. The more precise form for the top Yukawa coupling is 4> 2 / (M 2 + (p 2 )h,p v ^ v H, 
where h^ v ^ v n is a coupling of ipyipvH term, and <fi is the vev of (pi multiplied by 0i?/>Vv 
coupling. This is simply because the low energy third generation field is a linear combination 
of the form cos a ^3 — sin a ipv with the mixing angle sin a ~ 0/ \J Mf + <fr 2 . 

Therefore, in general, there is no gross contradiction to the fact that the top Yukawa 
coupling is order 1. However, in our case, if <fi/Mi becomes close to 1, the atmospheric 
mixing shifts from the maximal angle. Therefore, that needs to be addressed if the precise 
tri-bi-maximal mixing and h^ v ^ v n < 1 is demanded. The desired smallness of the effective 
/ and h! couplings however are more naturally obtained due to the presence of the Planck 
mass in the denominator. In order to make the /-coupling dominate over the h', we have to 
choose a small coupling for the H' Higgs field in Eq. (4). Similarly the A term in Eq. (9) is 
assumed to be small compared to the coefficient of the first matrix. 

Thus within these set of assumptions, this model is in good phenomenological agreement 
with observations. In a more complete theory, these assumptions need to be addressed. We 
however find it remarkable that despite these shortcomings, the model provides a very useful 
unification strategy of the diverse quark-lepton mixing patterns. 
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IV. VACUUM ALIGNMENT 



The major new point of this note is that we obtain the above fermion mass matrices from 
an S4 symmetry where the minimum configuration of the flavon fields used in our analyse of 
fermion mixings arise from superpotential minimization with very additional assumptions. 

We start our discussion by giving some simple examples and discussing the flavon align- 
ment as a prelude to the more realistic example. First thing to note is that 3i 3 is invariant 
under S4, but 32 3 is not. Denoting = (x,y, z), we see that in the first case, the singlet of 
3 = xyz. The superpotential for a 3i flavon field can therefore be written as 

W = ^m<p 2 - A0 3 = ^m(x 2 + y 2 + z 2 ) - Xxyz. (14) 

The solution of F-flat vacua (0 7^ 0) are 

= -{(1, 1, 1) or (1, -1, -1) or (-1, 1, -1) or (-1, -1, 1)}. (15) 

These aligned vacua can be identified to the vertex diagonal axes of the regular hexahedron. 
In fact £4 can be identified the permutation of the 4 axes of regular hexahedron. Once one 
of the axes is fixed, S'3-permutation is left. Therefore, the vacua break S4 down to S3. 

On the other hand, when 32 flavon is used (or the cubic term is forbidden by a discrete 
symmetry), quartic term involving the triplet is crucial for the F-flat vacua. The invariant 
quartic term 4 gives two linear combinations of the form x 4 + y 4 + z 4 and x 2 y 2 + y 2 z 2 + z 2 x 2 . 
This is because they have to be symmetric homogenous terms and invariant under the Klein's 
group, which is n rotation around the x, y, z axes. Thus, the superpotential term for 32 field 
is 

1 

W = -mtf - — (0 4 H - — (0 4 ) 2 (16) 

L ( 2 1 2 1 2\ / 4 1 4, 4\ / 2 2 1 22, 2 1\ 

= -far + y + z ) (x + y + z ) [x y + y z + z x ). 

The nontrivial F-flat vacua (0 7^ 0) are 



mM _ / mM -> / mM 

V k,W + 2 K ( 2 ) ' V ^5 + ^ C ' 

where a = (0,0, ±1), (0,±1,0), (±1,0,0), b = (±1,±1,±1), and c = (0,±1,±1), 
(±1,±1,0), (±1,0, ±1). We note that these vectors correspond to the axes of the regular 
hexahedron. The vacua break S4 down to Z4, Z3, and Z2, respectively. More importantly, 
the vacuum states in Eq. (7) used in the analysis of fermion masses in the previous section 
are a subset of the above vacua. 

Note that if we add a 4 term to the superpotential involving the 3i flavon field, a vacuum 
is possible, in addition to the original b vacua. However, c vacuum is absent. 
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Turning to the model at hand, due to non-trivial Z n charges for the flavon fields, the 
mass terms (bilinears) are not allowed. To solve this problem, we have included fields 
which then lead to Dirac type mass terms. The superpotential for a single flavon field is 
then given Wi = rriifafa + ^i^M- Denoting fa = (xi, y h Zi) and fa = (xi, y~i, zi). In terms of 
its component fields, we obtain 

Wi = mi(xiXi + y t yi + + -j^itftf + vlvl + z i z i) ( 18 ) 

K (2) _ _ _ K [3) 

Note that there are three kinds of invariant for fa 1 ^ 1 . Finding the F-flat solution of this 
superpotential is similar to the case in Eq.l fTTl) . It is easily verified that the F-flat vacua are 
proportional to a, b, and c similarly in Eq. (TTT1) . 
Several comments are now in order: 

• We note that the cubic terms in the flavon superpotential, such as 02^3 are forbid- 
den by our choice of Z n charge of Sj since their presence will spoil a vacuum alignment 
of 02 in Eq.©. 

• Secondly, note that the orthogonality of the vevs of <p2 and 03 is important to obtain 
the tri-bi-maximal mixing . One way to obtain it dynamically is to have a mixing 
term 0203 such that the coupling of the mixing term is much smaller than 0202 and 
0303 couplings. The invariant term 020| expressed in terms of components gives 
^2X32/22/3 + 2/2^/3^2^3 + ^2^3^2^ 3 , where 02 = (x 2 ,y2,z 2 ) and 3 = (x 3 ,y 3 ,z 3 ). The 
F-flatness condition implies that y 2 y 3 + z 2 z 3 = when x 2 = and x 3 ^ leading 
to the desired orthogonality of the alignments of (02) and (03). Note that with our 
Z n charge assignments, this can arise only in higher orders and its coefficients must 
therefore be small. The same situatrion happens also for the mixing terms of the form: 



>2<P2<P3<P3 



There are mixing terms between the different flavon fields in the quartic terms of 
the form Wij = -hfafa(f)j(ftj. When expressed in terms of the component fields 
x, y, z, they involve mixed terms like \(xiXiyji/j+yiyiZjZj +ZiZiXjXj) plus similar other 
mixed invariants. In the previous item, we just discussed the case when % = 2 and 
j — 3. As for the remaining terms of this type, they will in general induce small 
contributions proportional to A in the vevs in Eq.(7) where there are zeros. They 
will induce correction to the forms of our mass matrices. We will therefore need to 
assume that these A couplings to be small, so that their effect on our mass and mixing 
predictions will be small. 
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• Depending on the values of a and b, one could in principle get very high dimensional 
terms of the form s^s 2 4> 2 (p 2 4>3 (x,y are positive integers); however their contribution 
to the flavon potential is suppressed and we ignore these effects. 

We therefore conclude that all the desired vacua in the SO(10) model are present. Any 
possible corrections to them can be made small making it possible to take a first step towards 
building a unified model of flavor. 

V. CONCLUSION 

In summary, we have proposed a grand unified model for quark-lepton flavor starting 
above the GUT scale with an SO(10) theory with 5*4 x Z n discrete symmetry, S4 non-singlet 
flavon fields and two vector like pairs of 16 with mass above the GUT scale and SO(10) 
Higgs multiplets 10 and 126 fields that give mass to fermions. The 16 matter as well as the 
flavon fields transform as ^-family group triplets. The ground state of the flavon sector of 
the theory gives non-zero vevs to the flavon fields along specific directions due to the above 
discrete symmetries and when certain higher dimensional couplings between different flavon 
fields are assumed to be small. They fix the structure of the Yukawa couplings of 10 and 
126 fields at GUT scale after the vector-like fields decouple. This leads to specific mass 
textures for the quarks and leptons with only a few parameters and hence the predictions 
for quark lepton mass relations and mixing angles in both the quark and the lepton sector. 
In particular, the model leads to tri-bi-maximal form for the PMNS matrix in the leading 
order with corrections to this coming from charged lepton fields. Using this, we predict 
9 13 ~ 0.05. The quark mass hierarchies as well as quark mixings given by the model are in 
agreement with observations e.g. the model predicts at GUT scale correct mass ratios for 
mt,/m T and m s /m M as well as the Cabibbo angle V us without any adjustment of parameters. 
Some assumptions are needed to get the large top quark Yukawa coupling as well as relative 
strengths between the various flavon couplings. Clearly, our work begins a process which 
seems very promising and further work is needed to improve some of the assumptions used. 

Appendix : S4 group 

We briefly review the S4 group. The group S4 ~ D 2 x D 3 ~ (Z 2 x Z 2 ) x 
S 3 has irreducible reps li, 1 2 , 2, 3i and 3 2 as noted. To see the detailed 
properties, we use the (x, y, z) coordinate for the transformation law of the three- 
dimensional representations of S4. The group Z 2 x Z 2 is a Klein's group K = 
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{diag(l, 1, 1), diag(l, —1, —1), diag(— 1, 1, —1), diag(— 1, —1, 1)}, which corresponds to n ro- 
tation around the x,y,z axes. The group S 3 is a permutations of the three axes (x,y,z): 



/ n 1 n\ /nni \ / i n n 



S' = {diag(l J l,l) 



1 \ 

1 
1 

of 5*4 is given as S4 = {(k, s)\k G K,s G S*}. 

The 3i representation (column vector) transforms by the action of S4 as 



1 

1 

1 



V 



/ 



1 
1 
1 



1 

1 



1 

1 
1 



}. The element 



<p — > kscj), 

while 32 representation 0' transforms as 

cf)' — > (det s)kscf)' . 

The singlet I2 transforms as 



(19) 



(20) 



l a -> (dets)la, (21) 

and li is invariant under the action of S4. The reps li, I2, 2 are reps of S3 ~ D%, and the 
transformation law of 2 is rotation and reflection of the regular triangle. Doublet rep (u, v) 
transforms as 



VV 



->■ Uy B s U T b 



(22) 



For convenience, we list the Kronecker products of the triplets: 



(3i x 3i) 



2 ©3] 



(3i x 3|) a = 3 2 , 



3i x 3 2 = 1 2 © 2 ©3i ©3 2 . 
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